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ABSTRACT 

We consider a strong ultrarelativistic shock moving through a star whose 
envelope has a polytrope-like density profile. When the shock is close to the 
star's outer boundary, its behavior follows the self-similar solution given by Sari 
(2005) for implosions in planar geometry. Here we outline this solution and find 
the asymptotic solution as the shock reaches the star's edge. We then show that 
the motion after the shock breaks out of the star is described by a self-similar 
solution remarkably like the solution for the motion inside the star. In particular, 
the characteristic Lorentz factor, pressure, and density vary with time according 
to the same power laws both before and after the shock breaks out of the star. 
After emergence from the star, however, the self-similar solution's characteristic 
position corresponds to a point behind the leading edge of the flow rather than 
at the shock front, and the relevant range of values for the similarity variable 
changes. Our numerical integrations agree well with the analytic results both 
before and after the shock reaches the star's edge. 

Subject headings: hydrodynamics — shock waves — stars: general 



1. Introduction 

The surge of activity over the past decade or so in the fields of supernovae and of 
gamma-ray bursts and their afterglows has led to renewed investigation into the behavior of 
strong shocks. Much of the analytic work on strong shock propagation to date has focused 
on self-similar solutions to the hydrodynamic equations. In these solutions, the profiles of 
the hydrodynamic variables as functions of position have constant overall shapes whose time 
evolution consists simply of scalings in amplitude and position. As a result, self-similarity 
allows us to reduce the nominal system of two-dimensional partial differential hydrodynamic 
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equations to a system of one-dimensional ordinary differential equations. The existence of 
self-similar solutions thus enables a significant simplification of problems free of spatial scales 
in regions far from the initial conditions. The best-known such solutions are the pioneering 
Sedov- Taylor solutions for non-relativistic point explosions propagating into surroundings 
with power-law density profiles (Sedov 1946; von Neumann 1947; Taylor 1950). 

Self-similar solutions are traditionally divided into two categories (see, for example, 
Zel'dovich & Raizer (1967) for a detailed discussion). 'Type I' solutions are those in which 
the time evolution of the shock position and hydrodynamic variables follows from global 
conservation laws such as energy conservation. The Sedov- Taylor solutions are Type I; their 
ultrarelativistic analogues were found by Blandford & McKee (1976). By contrast, global 
conservation laws are useless in 'Type II' solutions, which are instead characterized by the 
requirement that the solution remain well-behaved at a singular point known as the 'sonic 
point'. If, for instance, the density of the surroundings falls off very quickly with distance, 
Type II solutions found by Waxman & Shvarts (1993) for non-relativistic spherical explosions 
hold instead of the Sedov- Taylor solutions and relativistic solutions found by Best & Sari 
(2000) hold instead of the Blandford-McKee solutions. 

Here we study the case of an ultrarelativistic shock wave moving outwards through a 
star whose envelope has a polytrope-like density profile. After the shock front reaches the 
outer edge of the star, an event we refer to as 'breakout', the shock front itself ceases to exist 
but the shocked fluid continues outward into the vacuum originally surrounding the star. 
We focus on the flow at times just before and just after breakout. As explained in §2, the 
shock evolution just inside the star's surface is identical to that expected for an imploding 
planar shock in a medium with a power-law density profile. Such a shock follows a Type 
II self-similar solution as discussed by Sari (2005) and Nakayama & Shigeyama (2005) and 
outlined briefly here. §3 describes the asymptotic solution as the shock front reaches the 
surface of the star, a singular point. In §4 we investigate the flow after breakout. We show 
that the self-similar solution for the evolution inside the star also describes the behavior 
outside the star except in that a different range of the similarity variable applies and in that 
the physical interpretation of the characteristic position changes. We show in §5 that the 
analytic results of §2, 3, 4 agree with our numerical integrations of the relativistic time- 
dependent hydrodynamic equations, and in §6 we summarize our findings. Throughout our 
discussion, we take the speed of light to be c = 1. 
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2. Shock propagation within the star 

Since we are interested in the shock after it has reached the envelope or the outermost 
layers of a star, we assume that the mass and distance lying between the shock front and the 
star's outer edge are much less than the mass and distance between the shock front and the 
star's center. In this region, we can take the star's gravity g to be constant and the geometry 
to be planar. We also assume that the stellar envelope has a polytrope-like equation of state, 
that is, p oc p q where p is the pressure, p is the mass density, and q is a constant. This type 
of equation of state occurs in various contexts including fully convective stellar envelopes, in 
which case q is the adiabatic index; radiative envelopes where the opacity has a power-law 
dependence on the density and temperature; and degenerate envelopes. 

Under these assumptions we can find the density profile from hydrostatic equilibrium 
and the equation of state as follows. Let x be the radial coordinate such that x = at the 
star's surface and x < inside the star. Then 

and with the boundary condition p = p = at the edge of the star, we have 

-^p*- 1 oc -gx (2) 

p oc (-xf/^ = (-x)~ k . (3) 

For convective and degenerate envelopes, q is between 4/3 and 5/3; for radiative envelopes 
with Kramers opacity, q = 30/17. These give k values between —1 and —3. 

With the power-law density profile p oc (— x)~ k , the evolution of an ultrarelativistic 
shock propagating through the envelope is given by a Type II converging planar self-similar 
solution to the hydrodynamic equations representing energy, momentum, and mass conser- 
vation, 

^b 2 (e + fp)}+-?-[ 1 ^(e + p)}=0 (4) 
|[7 2 /3(e+p)]+^[7 2 (/3 2 e + p)]=0 (5) 



with the ultrarelativistic equation of state 



P = 3 e • ( 7 ) 
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Here we will simply state the solution; for a detailed derivation see Sari (2005) or Nakayama 
& Shigeyama (2005). We assume the effect of the star's gravity on the shock propagation is 
negligible. Following Sari (2005), we let R(t) be the solution's characteristic position, which 
we choose to be the position of the shock front while the shock is within the star. We take 
t = at the time the shock reaches the star's surface (R — 0), and we take R < when 
t < 0. We take T(t), Pit), and N(t) to be respectively the characteristic Lorentz factor, 
pressure, and number density, and we define 



tT _ m tP tN _ m 

T~~2" ' -p-- m ~ k > aT-~7 



k 



Following Blandford & McKee (1976), we define the similarity variable as 



X = l + 2(m + l) 



R 



x 



R/T 2 



(8) 



(9) 



Note that for R < 0, x < R and the relevant range in x is — oo < x < 1 a s long as m > — 1. 
We define the hydrodynamic variables — the Lorentz factor 7, the pressure p, and the number 
density n — as follows: 

(10) 



1 \x ) t) = -T*(t)g{ X ) 



p(x,t) 
n(x,t) = N(t) 



P(t)f(x) 
Kx) 

7 l/2 (x) 



(11) 

(12) 



Here g, /, and h give the profiles of 7, p, and n; expressions for the dependence of m on k and 
for g, /, h as functions of x make up the entire self-similar solution. The above definitions 
and the ultrarelativistic hydrodynamic equations in planar geometry put the sonic point, 
the point separating fluid elements which can communicate with the shock front via sound 
waves from those which cannot, at gx = 4 — 2y / 3. Requiring that the solution pass smoothly 
through this point gives 

m = (3 - 2y/tj k 

(3-2\/3)fc 



9 = C g 



f = C 



9X 



2(2 + V3) 



3k-2kV3 + l 
-gX - 2k VS + 4 + 2VZ 



-2V3)k 



h = C h 



gX + 2*V3 - 4 - 2V3 



(2 v / 3-3)(2fc-l)fc 
(-l+kV3-V3) 



\9X 



2 -l+fcv^-v 7 ^ 



(13) 
(14) 
(15) 
(16) 
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The boundary conditions g(x = 1) = f(x — 1) — h(x = 1) = 1 which hold inside the star 
allow us to determine the constants of integration C g , Cf, Ch and write 



-gX - 2*V3 + 4 + 2^3 
-1-2*^3 + 4 + 2^ 



./ 



-9X 



2/c^3 + 4 + 2 v / 3 



h 



-1-2^ + 4 + 2^3 

(2V3-3)(2fc-l)fc 



(3-2x/3)ifc 



-(4-2V3)fc 



2^-4-2^3 



1 + 2/^3-4-2^ 



K 

[2 — gx] -i+^v^-v^ 



(17) 



(18) 



(19) 



3. Transition at breakout 

To know what happens to the shocked material after the shock front emerges from the 
star, we need the behavior of the shock just as the front reaches the surface — the 'initial 
conditions' for the evolution of the shock after breakout. Specifically, we are interested in 
the limiting behavior of each fluid element and in the asymptotic profiles of 7, p, and n as 
functions of x as t and R approach 0. 

The limiting behavior of a given fluid element may be found as follows. Due to the 
self-similarity, we know the time taken for 7, p, and n of a given fluid element to change 
significantly is the timescale on which R changes by an amount of order itself. Since R can 
change by this much only once between the time a given fluid element is shocked and the 
time the shock breaks out of the star, the limiting values of 7, p, and n for that fluid element 
should be larger only by a factor of order unity from their values when the fluid element was 
first shocked. 

We can also find the scalings of 7, p, and n with x at breakout via simple physical 
arguments. We denote by x , 70, po, n the position, Lorentz factor, pressure, and number 
density of a fluid element just after being shocked and by Xf, 7/, pf, rif those values at the 
time the shock breaks out. Since the shock accelerates to infinite Lorentz factors, and since, 
as we found above, the Lorentz factor of a given fluid element remains constant up to a 
numerical factor, this fluid element will lag behind the shock by Xf ~ Xq/^q at t = 0. Eq. 8 
gives T ~ t~ m / 2 , so we have 70 ~ (— x )~ m ^ 2 ; then 7/ ~ ((— Xf)^j)~ m ^ 2 or 

7/ ~ (_ X/ )--/2( m +i) _ (20) 
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Likewise, since P ~ t~ m ~ k and N ~ r m / 2 " fc , we have p ~ Xo m ~ fe and n ~ x m/2 fc ; then 

Pf ~ {{-x f )l))- m - k ~ (21) 
~ (-x/)- (m/2+fc)/(m+1) . (22) 

We can use the equations for the solution before breakout to perform equivalent calcu- 
lations of the limiting behavior of fluid elements and asymptotic profiles of 7, p, n. For the 
limiting behavior of a fluid element, we take the advective time derivative of gx and use the 
result to relate 7 and g to time for that fluid element. The advective derivative is given by 



D d „ d ■ d ■ d m + 1 _ . s d 

We apply this derivative to Eq. 17 to get 

D(g X ) (2 - g X ) (gx ~ 4 - 2^3 + 2y/3k) 



and integrate to get 
t/t 



9X 



Dlogt 



2 |(3+A/3)/(3fe-v / 3-3) 



(g X -*-2y/3) 



(23) 



(24) 



gx 



4-2 v / 3 + 2*V3 



1 - 4 - 2y/3 + 2kV3 



-3fc/(— \/3-3+3k) 



(25) 



where t is the time at which the fluid element is shocked, that is, when g = x = 1- When 
|<7X| ^ 1 — which becomes true everywhere behind the shock front as t — > — this simplifies 



to 

t/to 

and Eq. 17 simplifies to 



\9X\ 



1 - 4 - 2V3 + 2kV3 



3fe/(-v / 3-3+3fc) 



5 





-1-2*^3 + 4 + 2^. 



(3-2^3)*; 



(26) 



(27) 



We substitute Eq. 26 into Eq. 27 to get the limiting Lorentz factor of the fluid element as 
t -> 0: 

(3-3 v / 3)fe/(2(-\/3-3+3fc)) 



7 = 7o 



l-4-2v / 3 + 2fcv / 3 



(28) 



which is greater only by a numerical factor than the initial Lorentz factor 70 that the fluid 
element received right after being shocked. To relate the limiting p, n to p , n , we likewise 
take Eqs. 18, 19 in the limit \gx\ 3> 1 and use Eqs. 26, 27 with the results to get 



1 - 4 - 2^ + 2fc\/3 



(6-2^)fc/(- v / 3-3+3fe) 



(29) 
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n = no 



1 - 4 - 2x/3 + 2fcx73 



(4fe+fc v / 3-3-v / 3)(3-2\/3)fc 
2(fc- v /3-l-v / 3) 



(30) 



which again differ only by numerical factors from their values just after the fluid element is 
shocked. This is consistent with the behavior given above by simple physical considerations. 

For the analogous calculation of the asymptotic profiles of 7, p, and n, we cannot simply 
apply Eqs. 10, 11, 12: Eqs. 8, 9 require that x ~^ ~°° everywhere behind the shock and T, 
P, and iV diverge as t — > 0. Instead we take the t — > limit at a fixed position x. First we 
have 

X = I + 2(m + - x/R)T 2 ~2(m + l)(-x/R)T 2 . (31) 
With Eqs. 10, 27 this gives 



and 



7 





2 7 2 /r 2 


\ n (-Rr m - 


-l/2(l+m) - 


T 2 





9 



-A(m + l)(-x)^ 2 / R 
-1-2^ + 4 + 2^ 

4(m + 1) 



(32) 



-1-2^^ + 4 + 2x73. 



m/2(l+m) 



-m/2(l+m) 



(33) 



This is consistent with our qualitative discussion; the coefficient in the qualitative relation 
is a numerical factor times the constant (— R)~ m /T 2 . For the p and n profiles, we apply a 
similar analysis to the expressions for / and h in the limit t — > 0. 



p = P 



r 2 



m + fc 
1 + m 



4(m + 1) 



m-j-fc 
' 1 + m 



m+fc 

-X) 1+m 



n = iV 



r 2 



m/2+fc 

1 + m _ m/2+fc 

[4(m+ l)] __r +^ 



1-2^x73 + 4 + 2^/3 
_1 _ + 4 + 2x73] "*^ + - 1+W3 - V3 (- 



q/2+fc 



(34) 



•ro/2+fc 
1 + m 



(35) 



These results are likewise consistent with our qualitative discussion. 



4. Evolution after breakout 

4.1. Self-similar solution 

Since the breakout itself does not introduce new spatial scales into the flow, we expect 
the motion after breakout to remain self-similar. However, as the shock Lorentz factor 
diverges at t — 0, we cannot continue to associate the characteristic position, Lorentz factor, 
pressure, and number density with the values at the shock front after breakout. So we begin 
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by providing physical motivation for a different characteristic Lorentz factor and exploring 
the implications of this choice. 

We note that after breakout each fluid element expands and accelerates over time until 
the element's internal energy has been converted entirely into bulk motion. Given a rela- 
tivistic strong shock, the internal energy of a shocked fluid element in the frame moving with 
the fluid is comparable to the bulk kinetic energy of the fluid element. This implies that the 
fluid element's final bulk Lorentz factor should be much greater than the value of the shock 
Lorentz factor just after the fluid element was shocked. The timescale t x for the resulting 
expansion and acceleration is the time over which the fluid element's size and Lorentz factor 
change by a factor of order unity. For a fluid element located at —x and with Lorentz factor 
j x at t = 0, the time of breakout, this timescale is t x = x^l due to relativistic beaming. 
That every time t > is thus associated in a scale-independent way with a particular t x and 
7 suggests that we pick T(t = t x ) = ^ x to be the characteristic Lorentz factor. 

To see how V evolves with time, we use 7 oc (—^-"^(i+m) f rom gq_ 33 with the t x 
relation above to get T oc t"' m l 2 . For the characteristic pressure P and number density N, 
Eqs. 34, 35 likewise give P oc t~ m ~ k and N oc t~ m/2 ~ k . In other words, Eq. 8 holds after 
breakout with exactly the same k, m that apply inside the star. The characteristic position 
R is again the position which evolves according to the Lorentz factor V: R ~ 1 — l/2r 2 . 
Since the hydrodynamic equations still hold as well, Eqs. 9, 14, 15, 16 must remain valid 
when t > 0. 

To find the complete solution after breakout we need to specify the boundary conditions. 
We proceed by looking at the behavior of the similarity variables x-, 9-, f\ h. The relevant 
range in x depends on R, and while the relation between R and T is the same before and 
after breakout, R after breakout is not the position of the shock front. Instead, the front has 
infinite Lorentz factor and R lags further and further behind the shock with increasing time. 
A nice physical interpretation exists for R after breakout. R tracks the position corresponding 
to a fluid element which has expanded by a factor of order unity, so R marks the transition 
in position between fluid elements which have expanded and accelerated significantly since 
being shocked and fluid elements whose size and speed have remained roughly constant. 
Since it takes longer for fluid elements with smaller Lorentz factors to expand and accelerate 
significantly, R moves backward relative to the leading edge of the flow at x — t. Because 
R becomes positive after breakout, the range of possible x in the solution outside the star is 
x < t. Then x = at the 'front' x = t, and the relevant range in x m the solution outside 
is < x < 00 rather than —00 < x < 1. 

Far behind x = t, the profiles of 7, p, and n before breakout must coincide with the 
profiles after breakout. We know this because at a given time after breakout, sound waves 
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carrying the information that breakout occurred can only have traveled a finite distance 
behind the shock front; material further behind continues to flow as if the breakout had 
never occurred. Also, the two sets of profiles must coincide at t — 0, when everything is 
far behind the front. To phrase this requirement on the profiles in terms of the similarity 
variable, g(x — > — oo), f(x —> — oo)j and h(x - ► — oo) before breakout must coincide with 
g(x — > oo), fix — > oo), and h(x — > oo) after breakout. Then as x ^ oo after breakout, 
g, f, h — > and gx — > oo. In addition, the constants C 9 , Cf, Ch in Eqs. 14, 15, 16 must be 
the same for both the pre- and post-breakout solutions. In other words, the solutions before 
and after breakout, as specified by Eqs. 9, 14, 15, 16 and expressions for C g , Cf, Ch, are the 
same; only the relevant ranges in x an d the physical interpretations of the variables differ. 
So the expressions for g, /, h after breakout are 



9 = 



f 



gX + 2kV3-4:-2V3 
-1-2*^3 + 4 + 2^ 



9X 



+ 2*^3-4-2^ 



h 



-1-2^ + 4 + 2^3 
-gX - 2k V3 + 4 + 2^ 



(3-2 



-(4-2>/3)fc 



(2\/3-3)(2fc-l)fc 



K 

[gX — 2] -i+kV3-V3 



1 + 2^^-4-2^ 
The boundary conditions after breakout are given explicitly by g — f 



(36) 



(37) 



(38) 



1 and h = 

(5 + 4v^3 — 4^/3k) ^ ^+ k ^) a ^ ^ = 7 4_ z^y^ — 4\/3A;. A graphical comparison between 
the pre- and post-breakout 7 vs. position profiles is given in Figure 1 along with sample 
trajectories of fluid elements. 



4.2. Type I or Type II? 

While the flow before breakout follows a Type II self-similar solution, the solution 
describing the flow after breakout contains elements of Type I and Type II solutions. Unlike 
the Type II solution which applies before breakout, the post-breakout solution does not 
contain a sonic point. Differentiating Eq. 36 with respect to gx shows that the only local 
extremum of gx occurs at g = 00 or x — 0, where gx = 4 + 2^/3 — 2k\^3; since gx — > 00 
as x ~~ 00 ) 9X must attain its global minimum at x — 0- But then for k < neither the 
sonic point, gx = 4 — 2y/3, nor the other singular points, gx = 2 and gx = 4 + 2\/3, is 
included in the solution after breakout. A more physical argument for the exclusion of the 
sonic point from the post-breakout solution is that since each fluid element is accelerating 
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Fig. 1. — Profiles of 7 as a function of position (heavy lines) at seven different times marked 
on the figure and trajectories of three fluid elements in position-Lorentz factor space (thin 
lines). Fluid elements at the characteristic positions R are marked by open circles. We 
use x-t as the position coordinate to allow easy comparison of the profiles. The t — 
curve (heavy dotted line) is the asymptotic profile corresponding to the pure power law 
7 oc (_ a; )~"V 2 ( 1+m given in Eq. 33. The profiles with t < (heavy solid lines) are given 
by Eqs. 10, 17 and the profiles with t > (heavy dashed lines) are given by Eqs. 10, 36. 
When t < 0, the natural choices for R and T are respectively the location of the shock 
front and the Lorentz factor of the front. When t > 0, a fluid element at position R has 
accelerated by a factor of order unity and its Lorentz factor is of order V. So the positions 
R lie just above the 'knees' in the profiles, which separate fluid elements which have already 
expanded from those which have not. When \x — t\ ^> R/T 2 or, equivalently, \x\ ^00, all 
profiles approach the t — power law since at t — 0, \x\ — *• 00 everywhere behind the front. 
When \x — t\ <C R/T 2 , the t < profiles approach a constant (7 — > T/y/2) and the t > 
profiles approach 7 oc \x — (g oc x^ 1 from Eq. 36). Because every fluid element is always 
accelerating, the t < profiles always lie below the t — power law and the t > profiles are 
always above the t — power law. Trajectories of individual fluid elements before breakout 
are given by Eq. 25. After breakout, Eq. 25 still applies. The power laws relating t to gx 
stay the same after breakout since the equations for g before and after breakout are nearly 
identical; also, matching the pre- and post-breakout trajectories at t — gives the same \t \ 
in the evolution both before and after t — 0. 
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while T decreases with time, the fluid element moves forward relative to R and its x must 
decrease with time. Using Eq. 23 we see that D\j Dt < requires g\ > 2 > 4 — 2 a/3 for 
every fluid element. Then the entire post-breakout solution is causally connected as would 
be expected if it were Type I. 

Unlike Type I solutions, however, the solution after breakout contains infinite energy: 
it can be thought of as representing a flow into which a source at x = — oo feeds energy at 
a constant rate, sustaining the acceleration of fluid elements further and further behind the 
shock. As a result, global conservation laws do not apply just as would be expected in a 
Type II solution. So the post-breakout solution lies between the standard Type I and Type 
II solution categories. While this unusual situation implies that, in principle, the infinite 
energy contained in the solution can communicate with and affect the region near R, the 
regions of the solution containing this infinite energy lie arbitrarily far behind the front 
at x — t and therefore take arbitrarily long to communicate with the fluid near the front. 
Similarly, in any application of the post-breakout solution, the flow will be truncated at some 
position well behind R, potentially introducing a spatial scale into the problem. However, 
the solution is valid until information from the truncation region propagates to areas close 
to R. The further the truncation from R, the longer this will take. 



4.3. Relation to previous work 

The first analytic investigation of an ultrarelativistic planar shock wave was performed 
by Johnson & McKee (1971). The problem they consider is broadly similar to the one 
we discuss here, but our work differs in important respects from theirs. First, Johnson 
& McKee (1971) used the method of characteristics in their work: they analyzed the flow 
associated with the shock by tracing the paths of sound waves travelling through the fluid. 
Our analysis uses the self-similarity of the flow instead. So while some of their work can 
be applied to flows moving through fluids with arbitrary decreasing density profiles, their 
methods do not give profiles for the hydrodynamic variables as functions of a; at a given time. 
By contrast, our self-similar solutions require a power-law density profile inside the star but 
give explicit profiles for the hydrodynamic variables. Second, the methods used by Johnson 
& McKee (1971) require initial conditions consisting of a uniform stationary hot fluid about 
to expand into cold surroundings. In our scenario the hot expanding fluid is never uniform 
or stationary and always follows the self-similar profile specified by our solution. The self- 
similarity analysis tells us that the solution is Type II, at least before breakout; this implies 
that the asymptotic solution is independent of the initial engine. 

We can check that our asymptotic solution is consistent with the findings of Johnson & 
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McKee (1971) by looking at the Lorentz factors of individual fluid elements at very late times. 
While in our self-similar solution the fluid elements formally accelerate forever, each fluid 
element must in practice stop accelerating when all of its internal energy has been converted 
to bulk kinetic energy, or when p/n ~ 7///1 ~ 1. So we can estimate the final Lorentz factor 
of a given fluid element from Eqs. 36, 37, 38. By taking the advective time derivatives of 
7 and of ///t we can write differential equations for their time evolution following a single 
fluid element. These are 

L>7 
~Dt 

Dim 

Dt 




In the last steps we have taken the limit of late times when the accelerating fluid element 
approaches the shock front at % = 0. In this limit Eq. 36 implies g — > oo and gx — ► 
4 + 2\^3 — 2k\^3. Let 7 , /o, and h be the values of the functions in question just after our 
fluid element is shocked; then at late times 7 ^> 70 so (///i)/(/ //io) ~ l~ l ■ Integrating the 
above differential equations then gives 

/ , \ (V3-l)/2 

l=(f) ~7< 3 -^' 2 - T~7„ 1+ ^ ■ (42) 

7o \W 

This agrees with the results of Johnson & McKee (1971) for the final Lorentz factor of the 
fluid in a strong ultrarelativistic shock propagating into a cold medium with decreasing den- 
sity. The agreement provides additional support for our claim that the solution outside the 
star behaves like the solution describing a standard planar shock up to the initial conditions 
and the interpretation of the characteristic values R, T, P, N. Note that the differences be- 
tween the initial conditions used in their work and in ours are unimportant to the scaling law 
relating the final and initial Lorentz factors of a given fluid element. This result agrees with 
the findings of Tan et al. (2001) concerning the scaling law: partly because of uncertainty 
over the different initial conditions, they used numerical simulations to check the 7 ~ 7o + ^ 
result. 

Recently, Nakayama & Shigeyama (2005) also investigated the problem of an ultrarela- 
tivistic planar shock. While the self-similar solution they give for the flow before breakout 
is identical to the one in Sari (2005) and outlined here, they do not give analytic results for 
or a physical interpretation of the self-similar solution after breakout. 
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5. Comparison with numerical integrations 

To verify our results numerically, we integrated the time-dependent relativistic hydro- 
dynamic equations using a one-dimensional code. Figure 2 shows curves for 7 as a function 
of position at a single time before breakout, while Figure 4 shows the time evolution of T, 
P, and N before breakout. The numerical and analytic results are in excellent agreement. 
Figures 3 and 5 respectively show the 7 vs. x profile and time evolution of T, P, and N after 
breakout; the agreement between numerical and analytic results here confirms the choice of 
scale R(t) after breakout that we discussed in § 4.1. 

6. Summary 

We have shown that, given an ultrarelativistic shock propagating into a planar polytropic 
envelope, the flow upon the shock's emergence from the envelope into vacuum follows a self- 
similar solution strikingly similar to the self-similar solution describing the flow while the 
shock remains within the envelope. Both self-similar solutions obey the same relations with 
regard to the time-evolution of the characteristic quantities R, T, P, N and with regard to 
the similarity variables \i fiS /> h. The pre- and post-breakout solutions differ only in that 
the applicable ranges in x an d the physical interpretations of the characteristic quantities 
differ. As a result of these differences, the behavior of the flow after breakout lies somewhere 
between the traditional Type I and Type II classes of self-similar solutions; before breakout 
a Type II solution applies. To arrive at these results we have looked in detail at the behavior 
when the shock reaches the outer edge of the envelope. 

We have discussed these results in the context of an application — the motion of a shock 
wave through a polytropic envelope near the surface of a star, the shock's emergence from 
the surface, and the subsequent flow into vacuum. This situation may be related to the 
explosions believed to cause gamma-ray bursts and supernovae (see, for example, Tan et al. 
2001) and should be especially relevant in very optically thick media such as neutron stars. 

This research was partially funded by a NASA ATP grant. RS is a Packard Fellow and 
an Alfred P. Sloan Research Fellow. 
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Fig. 2. — Lorentz factor 7 as a function of position x shortly before the shock breaks out of 
the star. The density profile has power-law index k = —1.5. The analytic profile taken from 
the self-similar solution (solid line) agrees well with the numerical profile (crosses). 
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Fig. 4. — Evolution of T (top panel), P (middle panel), and N (bottom panel) with R while 
the shock is still inside the star. The density profile has power-law index k = —1.5. The 
evolution of T, P, N with R is equivalent to time evolution when T ^> 1. Crosses represent 
numerical data; solid lines are the best-fit lines to the data. That the data are well fit by lines 
implies that T, P, and iV do indeed evolve as power laws; that the numerical and analytic 
slopes agree confirms that the evolution is as expected. 
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Fig. 5. — Same as Figure 4 but for times after the shock emerges from the star. T, P, N, 
and R were deduced from the numerical data by finding at each time the position where 
7~ 1 p/n, the ratio of the thermal and the bulk kinetic energies in the frame of the fluid, fell 
just below the constant value we expect at the time of breakout. 
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